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A  neighboring  extremal  guidance  law  for  systems  with  piecewise  linear 
control  using  multiple  optimization  along  the  suboptimal  trajectory  is  devel¬ 
oped.  The  guidance  law  computes  control  gains  that  relate  the  difference  in  the 
perturbed  and  nominal  trajectory  with  a  control  update.  Previous  research  in 
neighboring  extremal  control  has  concentrated  on  calculating  the  control  gains 
using  a  single  optimization  from  the  initial  conditions  to  the  final  constraint 
manifold.  The  purpose  of  this  study  is  to  develop  a  more  optimal  method 
of  computing  the  control  gains  at  each  node  through  multiple  optimization. 

o\9^2-31212 

IHHHIRIIII^ 


92  12  09  038 


Multiple  optimization  reoptimizes  the  trajectory  from  each  control  node  on 
the  suboptimal  trajectory.  The  suboptimal  control  vectors  for  each  subtra¬ 
jectory  are  used  to  compute  a  more  optimal  set  of  control  gains.  An  example 
optimization  problem,  the  lunar  launch  problem,  is  used  to  compare  the  perfor¬ 
mance  of  both  the  single  optimization  and  multiple  optimization  control  laws. 
Overall,  the  multiple  optimization  gains  returned  results  comparable  to  those 
achieved  with  single  optimization.  For  perturbations  of  ±5%  in  either  grav¬ 
ity  or  thrust  acceleration,  endpoint  conditions  for  multiple  optimization  gains 
never  exceeded  more  than  50  feet  in  position  or  4.3  ft/sec  in  vertical  velocity 
from  the  desired  boundary  conditions.  Direct  comparisons  of  performance  for 
specific  model  perturbations  result  in  mixed  conclusions.  For  thrust  pertur¬ 
bations,  single  optimization  gains  deleiver  smaller  errors  in  meeting  end  point 
conditions,  while  for  gravity  perturbations,  multiple  gains  give  better  results. 
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Chapter  1 


Introduction 


An  inherent  problem  in  optimal  and  suboptimal  control  is  that  the 
optimal  trajectory  is  not  usually  realizable.  Unpredicted  model  errors  and/or 
perturbations  in  the  initial  conditions  can  result  in  deviations  from  the  optimal 
trajectory.  Recomputing  the  revised  optimal  control  from  the  point  of  pertur¬ 
bation  to  the  final  constraint  manifold  is  not  always  possible  due  to  computer 
processing  limitations.  As  a  result,  the  precomputed  optimal  control  does  not 
produce  the  desired  trajectory. 

In  light  of  these  deviations,  one  approach  to  optimal  control  is  to  pro¬ 
vide  guidance  in  the  neighborhood  of  the  optimal  trajectory.  Assuming  model 
errors  and  perturbations  are  small  enough,  the  neighboring  extremal  guidance 
should  keep  the  perturbed  trajectory  in  the  neighborhood  of  the  extremal,  or 
optimal,  path. 

The  purpose  of  this  research  is  to  develop  a  multiple  optimization 
neighboring  extremal  guidance  law  for  systems  using  piecewise  linear  control 
and  apply  it  to  a  representative  problem.  The  guidance  law  provides  updates  to 
the  suboptimal  control  vector  based  on  errors  in  the  perturbed  versus  the  nom¬ 
inal  trajectory  and  precomputed  neighboring  extremal  gains.  The  neighboring 
extremal  control  gains  are  found  by  finding  the  change  in  control  required  to 
minimize  the  change  in  the  performance  index  subject  to  the  trajectory  reach- 
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ing  the  final  constraints. 

Two  methods  of  computing  these  gains  are  presented  and  discussed. 
In  the  first  method  the  gains  are  computed  along  the  optimal  trajectory  at  node 
points  corresponding  to  the  piecewise  linear  control  nodes.  The  second  method 
reoptimizes  from  each  node  along  the  optimal  trajectory  to  the  final  constraint 
manifold.  The  neighboring  extremal  gains  are  then  computed  by  integrating 
back  from  the  final  constraint  manifold  along  the  reoptimized  trajectory  to 
each  node. 

To  test  each  approach,  a  simple  example,  the  lunar  launch  problem,  is 
used.  The  lunar  launch  problem  involves  placing  a  vehicle  in  a  low  moon  orbit, 
in  minimum  time,  subject  to  constraints  on  the  final  conditions.  Perturbations 
in  vehicle  acceleration  and  lunar  gravity  are  used  to  test  the  ability  of  the 
control  law  to  correct  for  “model  errors”  and  reach  the  desired  final  conditions. 


Chapter  2 


Suboptimal  Control  Problem 


The  fixed  final  time  problem  is  to  find  the  control,  u(t),  that  minimizes 
the  performance  index 


J  -  <j>(xj), 


(2.1) 


with  system  dynamics 

x  =  f(t,x,u), 

and  initial  conditions 

to  —  to.  i  xo  =  *o., 
subject  to  boundary  constraints 


(2.2) 

(2.3) 


0(x/)  =  0. 


(2.4) 


The  solution  to  this  system  of  equations  produces  the  optimal  control 
history,  tt(f).  However,  often  the  system  of  equations  cannot  be  solved  in  a 
closed  form  or  may  be  too  complicated  to  solve  numerically.  By  parametriz¬ 
ing  the  control,  a  nonlinear  programming  code  is  employed  to  solve  for  the 
suboptimal  control. 

Consider  a  suboptimal  control  that  is  of  the  class  of  piecewise  linear 
functions.  For  this  study  only  a  single  control  is  considered;  however,  the  results 
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can  be  generalized  to  include  multiple  control  systems.  Hence,  the  control  is 
expressed  as  a  parameter  vector 

aT  =  [an,ai2>,-'»air]5  (2.5) 

where  r  is  the  number  of  control  points  or  nodes  along  the  trajectory.  By 
normalizing  the  variable  of  integration  so  the  integration  limits  are  0  to  1,  the 
control  nodes  become  fixed  as  fractional  representations.  If  a  new  variable  of 
integration,  r,  is  defined  as 

t  =  f ,  (2.6) 

*/ 


the  system  dynamics  become 


where 


x'  =  —  =  tjf(tfT,x,a)  =  g(r,tf,x,a), 


t0  =  0  ,  Tf  =  1. 


Since  the  states  are  a  function  of  the  parametrized  control,  the  fixed 
final  time,  suboptimal  control  problem  is  to  find  the  piecewise  continuous  con¬ 
trol,  a,  that  minimizes  the  performance  index 

J  =  «*,(<>)),  (2.9) 


with  boundary  constraints 


^(x/(a))  =  0, 


subject  to  system  dynamics 


(2-10) 


(2.11) 
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where 

t0  =  0  ,  Tf  =  1.  (2.12) 

The  suboptimal  control  vector,  a,  is  obtained  using  the  nonlinear  pro¬ 
gramming  code  VF02AD.  Equipped  with  the  suboptimal  control  and  trajectory, 
the  neighboring  extremal  guidance  law  can  be  developed. 


Chapter  3 


Neighboring  Extremal  Control 


The  suboptimal  control  solution  provides  a  control  history  and  a  nom¬ 
inal  trajectory  that  minimizes  the  performance  index.  As  long  as  the  vehicle 
remains  on  the  nominal  trajectory,  the  control  law  guides  the  vehicle  to  the 
desired  final  conditions  in  the  minimum  time.  Realistically,  however,  the  nom¬ 
inal  trajectory  is  rarely  achieved.  Perturbations  in  initial  conditions  or  model 
errors  produce  deviations  from  this  trajectory.  The  goal  of  neighboring  ex¬ 
tremal  control  is  to  develop  a  guidance  law  for  the  perturbed  trajectory  in  the 
neighborhood  of  the  nominal  trajectory.  The  guidance  law  provides  control 
corrections  that  ultimately  reduce  the  change  in  the  performance  index  of  the 
perturbed  trajectory  from  that  of  the  nominal  trajectory. 

3.1  Neighboring  Control  Derivation 

The  basis  of  developing  a  neighboring  extremal  control  law  is  to  find 
the  change  in  control  required  to  minimize  the  increase  in  the  performance 
index  from  the  point  of  the  perturbation  to  the  final  constraint  manifold.  The 
derivation  of  this  control  law  was  originally  developed  by  Hull  and  Helfrich  [1]. 
A  summary  of  the  basic  derivation  is  included  here  for  completeness. 

Since  xt  —  xj(a)  from  Eq.  (2.11),  the  augmented  performance  index 
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for  the  suboptimal  control  problem  is  expressed  as 

J'  =  +  vTrl>(xf)  =  G(a,  i/),  (3.1) 

where  <f>  is  the  performance  index,  0  is  the  vector  of  final  constraints  and  v  is 
the  vector  of  Lagrange  multipliers.  The  change  in  the  performance  index  for  a 
perturbation  in  the  control  vector,  a,  to  second  order  is 

A  J'  =  GaSa  +  1 6aT(%a6a  (3.2) 

Since  the  first  variation  vanishes  on  the  extremal  path,  the  change  in  the  per¬ 
formance  index  is  approximated  as 

AJ'  =  ±6aTGla6a.  (3.3) 


In  Eq.  (3.3),  the  6a’ s  are  not  independent  since  the  final  conditions 
must  be  satisfied.  Hence,  it  is  necessary  to  determine  the  relationship  between 
6a  and  6x o  such  that  6rp  =  \})Xj6xj  =  0.  Eq.  (2.11)  relates  the  state  and  the 
control  through  a  set  of  differential  equations.  Taking  the  variation  of  this 


equation  results  in 


with  boundary  conditions 


A  solution  of  the  form 


6x'  =  gx6x  +  ga6a, 

(3.4) 

fH 

II 

c* 

o 

II 

»? 

(3.5) 

O 

II 

t 

• 

(3.6) 

ipXf6xj  =  0. 

(3.7) 

6x  =  $6x/  + 

(3.8) 

8 


is  assumed,  with  boundary  conditions 

•-T 

II 

•e> 

(3.9) 

*/  =  o, 

(3.10) 

and  substituted  into  Eq.  (3.4).  This  results  in 

($'  -  gx$)6xf  +  (<&'  -  gxty  -  ga)6a  =  0. 

(3.H) 

To  satisfy  Eq.  (3.11)  the  coefficients  of  6xj  and  6a  are  chosen  to  equal  zero. 
This  produces  the  differential  equations  governing  $  and 

=  (3.12) 

and 

*'  =  **  +  *  (3-13) 

with  final  conditions 

*,  =  /,  (3.14) 

*/  =  0.  (3.15) 

These  differential  equations  are  integrated  backward  from  tj  to  find  $  and 
at  any  node.  Multiplying  Eq.  (3.8)  by  $-1,  solving  for  6xj,  and  substituting 
into  Eq.  (3.7)  gives  the  desired  relation: 

-  V’r/^o1^  =  0.  (3.16) 

This  boundary  constraint  is  combined  with  the  performance  index  in 
Eq.  (3.3)  to  form  a  new  augmented  performance  index 

J'  =  +  SvT(f.,tZ,Vofo  -  4>x/*o,{xo)- 


(3.17) 
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Minimizing  J'  minimizes  the  change  in  the  perfonnance  index  due  to  a  per¬ 
turbation  in  the  control,  subject  to  the  boundary  constraint  that  relates  the 
control  and  the  state  vector.  The  matrix  Gaa  is  not  necessarily  positive  definite 
since  6a  has  both  dependent  and  independent  variations.  However,  the  bound¬ 
ary  constraint  is  used  to  eliminate  the  dependent  variations  so  the  resulting 
matrix  is  composed  of  only  independent  variations,  and  the  performance  index 
is  in  quadratic  form. 

The  first  variation  of  this  performance  index  is 

6J'  =  (6aTGaa  +  6vTrl>X/$ol*0)6(6a).  (3.18) 

To  satisfy  the  requirement  that  the  first  variation  vanish,  the  coefficient  of  6(6a) 
must  be  zero.  Therefore, 

6aTGaa  +  6vTtl>Xl$  =  0.  (3.19) 

Solving  for  6a, 

6a  =  -G2*l*ZT+xt6v  (3-20) 

substituting  into  the  constraint,  Eq.  (3.16),  and  solving  for  6v  gives 

6v  =  -(xJ>X/$o1V0GA\VZ$oTrl;Zj)-1rJ>x/$ol6xo  .  (3.21) 

Hence,  Eq.  (3.20)  leads  to 

6a  —  K06x0,  (3.22) 

where 

K»  =  (3.23) 

is  the  gain  matrix  at  To.  Eq.  (3.22)  relates  the  change  in  the  control  required  to 

minimize  the  change  in  the  performance  index  for  a  perturbation  in  the  state 
at  Xq  and  still  satisfy  the  final  conditions. 
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3.2  Single  Trajectory  Optimization 

For  continuous  systems  the  neighboring  extremal  guidance  law  can 
be  applied  throughout  the  trajectory.  By  the  optimality  condition  [2],  the  tra¬ 
jectory  from  each  point  along  the  optimal  path  to  the  final  constraint  manifold 
can  be  thought  of  as  the  optimal  trajectory  from  that  point.  As  a  result,  the 
gain  equation  for  the  continuous  optimal  control  problem  is  valid  throughout 
the  trajectory. 

For  the  suboptimal  problem,  however,  the  optimality  condition  does 
not  apply.  If  the  suboptimal  trajectory  has  r  control  nodes,  the  suboptimal 
trajectory  from  node  1  is  an  r-node  optimal  trajectory.  FYom  the  second  node 
though,  the  suboptimal  trajectory  is  an  (r  —  l)-node  suboptimal  trajectory. 
Furthermore,  in  order  to  satisfy  a  suboptimal  control  problem  with  m  bound¬ 
ary  constraints,  at  least  m  control  nodes  must  be  available.  For  a  parametrized 
control  with  r  control  nodes,  m  nodes  are  necessary  to  solve  the  boundary  con¬ 
ditions  while  the  remaining  r  —  m  nodes  (m  <  r)  provide  the  optimization.  As 
a  result,  neighboring  extremal  gains  can  only  be  computed  along  the  subopti¬ 
mal  trajectory  at  the  first  r  —  ra  +  1  node  points.  The  remaining  m  —  1  node 
points  do  not  have  enough  control  nodes  to  satisfy  the  boundary  constraints. 

The  single  trajectory  optimization  method  computes  the  neighboring 
extremal  gains  by  evaluating  the  gain  equation  (Eq.  3.23)  with  the  values  of 
$  and  $  along  the  suboptimal  trajectory.  The  differential  equations,  $'  and 
are  integrated  back  from  the  final  constraint  manifold  to  each  node  on  the 
suboptimal  trajectory  using  the  suboptimal  control  vector.  Gain  matrices  are 
computed  at  each  node  point  up  to  node  r  —  m  +  1.  Since  the  calculation  of 
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Kp  is  based  on  the  nominal  trajectory,  the  gains  can  be  precomputed  at  node 
points  and  stored. 

Although  the  gain  matrix  at  each  node  is  am  rxn  matrix,  where  r 
is  the  number  of  nodes  and  n  are  the  number  of  states,  only  gains  associated 
with  the  current  node  are  retained.  Each  gain  matrix  relates  the  change  in 
control  along  the  suboptimad  trajectory  to  am  error  in  the  state  vector  at  a 
node  point  adong  the  trajectory.  Obviously,  control  changes  at  nodes  that  have 
been  “passed”  do  not  affect  the  current  control  amd  are  eliminated.  Similarly, 
gains  related  to  nodes  “in  front”  of  the  current  node  are  also  eliminated.  As  a 
result,  the  matrix  reduces  to  a  1  x  n  vector. 

Since  each  node  adong  the  trajectory  corresponds  to  a  particular  r,  the 
gain  vector  at  auiy  rp  is  found  by  lineaudy  interpolating  adjacent  gain  vectors. 
As  the  last  gain  vector  is  at  node  r  —  m  -f  1,  gain  vectors  beyond  this  node  are 
linearly  extrapolated  based  on  the  last  two  gain  vectors. 
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3.3  Multiple  Trajectory  Optimization 

The  focus  of  this  research  has  turned  to  developing  another  method 
of  computing  the  control  gains  along  the  suboptimal  trajectory.  This  approach 
uses  reoptimized  r-node  control  vectors  from  each  node  along  the  suboptimal 
trajectory  to  compute  the  gain  matrices.  Fig.  3.1  depicts  a  five-node  example 
of  a  trajectory  and  subtrajectory.  To  prevent  confusion,  the  trajectory  from 
any  node  to  the  final  constraint  manifold,  except  for  the  first  node,  is  called 
a  suboptimal  subtrajectory.  The  suboptimal  subtrajectory  from  the  first  node 
to  the  final  constraint  manifold  is  also  the  suboptimal  trajectory. 

Each  subtrajectory  control  vector  is  slightly  different  from  the  sub¬ 
optimal  trajectory.  As  expected  the  suboptimal  trajectory  from  node  1  is  the 
same  as  in  the  single  optimization  method.  However,  from  the  second  node, 
multiple  optimization  provides  an  r-node  subtrajectory  to  the  final  constraint 
manifold,  while  the  single  optimization  method  only  has  r  —  1  nodes  from  the 
second  node  to  the  final  constraint  manifold. 

In  the  multiple  optimization  approach,  the  differential  equations  for  $ 
and  are  integrated  back  to  each  node  using  the  optimized  control  vector  from 
that  node.  Only  the  gain  matrices  at  nodes  along  the  suboptimal  trajectory  are 
stored.  In  the  case  of  node  1,  the  suboptimal  control  vector  is  used  to  integrate 
$  and  ¥  back  from  the  final  constraint  manifold  to  node  1.  For  a  subsequent 
node,  the  optimal  control  vector  along  the  subtrajectory  is  used  to  integrate  $ 
and  ^  back  to  the  node.  This  process  is  repeated  for  the  first  r  —  m  +  1  nodes 
along  the  suboptimal  trajectory. 

The  gain  matrix  at  each  node  is  reduced  to  a  lxn  vector  by  eliminat- 
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Figure  3.1:  Optimal  trajectory /subtrajectory 

ing  gains  that  do  not  affect  the  control  along  the  suboptimal  trajectory.  Since 
the  neighboring  extremal  path  is  based  on  the  original  suboptimal  trajectory, 
only  those  gains  that  affect  the  control  along  the  suboptimal  trajectory  need 
to  be  retained.  Only  the  first  row  of  the  gain  matrix  at  each  node  relates  to 
the  control  along  the  trajectory. 

The  hypothesis  behind  optimizing  from  each  node  point  is  that  the 
resulting  gains  are  based  on  the  most  optimal  r  node  trajectory  from  any  node 
along  the  suboptimal  trajectory.  In  the  single  optimization  method  the  gain 
matrices  are  based  on  the  original  suboptimal  trajectory  from  node  1,  and 
the  control  vector  updates  are  also  along  this  trajectory.  Since  the  multiple 
optimization  approach  bases  its  gains  on  the  optimal  r  node  subtrajectory 
from  each  node  point,  the  gain  vector  should  produce  a  control  update  that  is 
more  aligned  with  the  most  optimum  path  to  the  final  constraint  manifold  but 
still  in  the  neighborhood  of  the  extremal. 


Chapter  4 


Lunar  Launch  Problem 

4.1  Optimal  Control  Problem 

The  test  problem  for  this  research  involves  injecting  a  vehicle  into 
lunar  orbit  in  minimal  time  (refer  to  Fig.  4.1).  The  nominal  solution  to  this 
free  final  time  problem  produces  a  control  history  (thrust  angle  versus  time) 
that  minimizes  the  performance  index 

J  =  tf.  (4.1) 

The  differential  equations  governing  the  trajectory  are 

y  =  t>  (4.2) 

u  =  acosff  (4.3) 

v  =  a  sin  0  —  g,  (4.4) 

where  a  is  the  vehicle  acceleration  and  g  is  the  acceleration  of  gravity.  The 
differential  equations  for  the  downrange  component  are  not  included  since  they 
are  not  constrained.  The  initial  conditions  are  given  by 

t0  =  0,  yo  =  0,  uo  =  0,  u0  =  0  ,  (4.5) 

and  the  final  conditions  by 

y,  =  50,000  ft  (4.6) 
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Figure  4.1:  The  lunar  launch  problem 

uj  =  5,444  ft/ sec 
Vf  =  0  ft/ sec  . 


These  final  conditions  are  the  desired  final  states  to  achieve  lunar  orbit.  The 
quantities  u  and  v  denote  the  horizontal  and  vertical  components  of  the  velocity, 
respectively.  The  final  boundary  constraints  are  combined  into  a  vector,  rp, 


where 


50000  1 

=  ^4-1  =0. 


For  simplicity  several  assumptions  are  made.  First,  the  thrust  vector 
and  all  motion  is  in  the  x-y  plane.  Second,  it  is  assumed  that  the  moon  is 
flat  over  the  course  of  the  trajectory  so  the  lunar  gravity  vector,  g,  remains 
constant  with  j<?|  =  5.32  ft/sec2.  Finally,  the  ratio  of  vehicle  mass  and  thrust 
is  assumed  constant  throughout  the  trajectory  so  thrust  acceleration,  a,  is 
constant.  Based  on  previous  research  into  the  lunar  launch  problem  [1],  a 
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nominal  value  of  a  =  20.8  ft/sec2  is  assumed. 

In  its  current  form  as  a  free  final  time  problem,  the  upper  integration 
limit  of  the  equations  of  motion  is  not  fixed.  By  selecting  a  different  indepen¬ 
dent  variable,  the  integration  limit  can  be  fixed,  and  the  problem  handled  more 
conveniently. 

For  the  lunar  launch  problem,  the  control  angle,  0,  stays  between  0 
and  90  degrees  during  the  ascent.  As  a  result,  u  is  monotonic  and  it  is  used  as 
the  variable  of  integration. 

With  u  as  the  variable  of  integration,  several  benefits  are  realized. 
Since  u/  is  the  integration  limit,  the  vehicle  is  certain  to  meet  the  required 
horizontal  velocity  to  achieve  orbit.  In  addition,  incorporating  u  into  the  limits 
of  integration  reduces  the  number  of  boundary  constraints  since  uj  will,  by 
definition,  be  satisfied. 

The  state  vector  is  rewritten  as 

xT  =  [  t  y  v  ] .  (4-10) 

To  normalize  by  uj  the  variable  u  is  introduced  where 

U  —  Uo 

u  = - . 

Uf-Uo 

For  the  lunar  launch  problem,  «/  =  5444  ft/sec  and  u<j  is  the  initial  horizontal 
velocity  of  the  vehicle.  For  simulations  and  optimization  of  the  nominal  trajec¬ 
tory  from  u  =  0,  «o  equals  zero.  For  the  single  trajectory  optimization  uo  =  0 
during  gain  calculations,  while  in  the  multiple  optimization  scheme,  uo  equals 
the  horizontal  velocity  at  each  control  node. 


(4.11) 
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The  normalized  differential  equations  are 


dx 

du 


dx  dt  du 
dt  du  du 


g(x,u,a). 


Therefore,  the  equations  of  motion  are  rewritten  as 


t,  _  («/  ~  “o) 
a  cos  9 

s'  = 

a  cos  9 

v,  _  (a  sin  9  -  g)(uj  -  u0) 
a  cos  9 

where 

Xq  =  [  0  0  0  ] , 

y/  =  50000, 


Vj  =  0 


and  the  limits  of  integration  are 

«o  =  0  ,  S/sl. 


(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 


The  suboptimal  control  vector  is  computed  using  the  nonlinear  pro¬ 
gramming  code  VF02AD.  VF02AD  is  an  extension  of  a  variable  metric  method 
for  unconstrained  optimization  applied  to  the  constrained  case.  It  searches  for 
an  optimal  solution  until  reaching  a  preset  convergence  criteria.  Convergence 
occurs  when  the  change  in  the  performance  index  between  iterations  is  less 
than  the  preset  criteria.  For  all  optimizations  a  convergence  criteria  of  10~9  is 
used. 
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For  the  single  trajectory  optimization,  only  the  suboptimal  control 
from  tio  =  0  to  Uf  =  5444  ft /sec  is  computed.  The  suboptimal  control  consists 
of  nine  nodes  with  an  interval  of  680.5  ft/sec  between  node  points.  Nine  nodes 
are  selected  for  convenience,  resulting  in  control  nodes  every  0.125u. 

For  the  multiple  optimization  method,  the  nine-node  suboptimal  tra¬ 
jectory  from  u  =  0  is  computed  and  used  as  a  basis  for  the  nominal,  unper¬ 
turbed  control.  Then  the  nine-node  suboptimal  subtrajectory  from  each  node 
on  the  suboptimal  trajectory  to  the  final  constraint  manifold  is  calculated. 
The  control  vector  for  each  subtrajectory  gives  the  suboptimal  nine-node  con¬ 
trol  from  tio  =  680. 5(N  —  1)  ft/sec  to  ti/  =  5444  ft/sec,  where  N  is  a  node  on 
the  suboptimal  trajectory.  This  results  in  eight  suboptimal  control  vectors  for 
the  multiple  optimization  approach. 

4.2  Computation  of  the  Neighboring  Extremal  Gains 

To  compute  the  neighboring  extremal  gains,  expressions  for  Gaa,  tpX/, 
and  $  are  found. 

The  derivative,  Gaa ,  is  calculated  numerically  using  central  differ¬ 
ences.  Gaa  represents  the  second  derivative  effect  of  a  change  in  the  parametrized 
control  vector  on  the  performance  index,  G.  Each  element  of  a  is  perturbed 
and  the  effect  on  the  performance  index  is  computed. 

tl>Xf  is  the  partial  derivative  of  the  constraints  with  respect  to  the 
state  vector.  rj>Xf  is  a  constant  matrix  with 

tx,  = 


(4.20) 
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The  differential  equations  for  $  and  *P  are  integrated  backwards  from 
tt /.  The  differential  equation  governing  $  is 

d$  .  _ 


=  *'  =  gt*, 


(4.21) 


where 


,g>|  H 
|C£> 

II 

'  0  0 

0  ■ 

=  00 

U/-U0 
a  cot  8 

(4.22) 

.0  0 

0 

p,  =  / 

(4.23) 

The  matrix  differential  equation  for  P  is 

P'  =  s*P  +  9a 

where 


(4.24) 


P/  =  0.  (4.25) 

In  tb;s  problem,  P  is  a  3x9  matrix  and  ga  is  the  partial  derivative  of  the  per¬ 
formance  index  with  respect  to  the  parameter  vector  a.  Since  aT  =  [0 i,  •  •  • ,  0r]» 
ga  can  be  written  as 


(4.26) 


(«/  —  uo)  sin  6  si  "Si  ,3 

Ot  COS2  9  (a-jiinl)  gg  (q-gain?)  99  (q-g  ting)  98 

•  am 0  un(  d& 2  sin 0  dfl#  • 

The  control,  0,  is  linearly  interpolated  between  adjacent  nodes,  so  that 

0-  —  0 

0(u)  =  0,_!  +  — - =^(0  -  u.-x)  ,  for  Uj_i  <  u  <  u,, 

U,  —  Uj_i 


(4.27) 


0(tt)  =  0;  +  — :r(ti  —  «i)  ,  for  u;  <  u  <  u,+i.  (4.28) 

u.+i  -  «. 
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The  partial  derivatives  with  respect  to  respectively,  are 

dO  _  u  — 

dOi  Hi  -  u,_i  ’ 


(4.29) 


and 


dO  _  €t  —  tx< 

dOi  tti+i  -  Hi  ’ 


(4.30) 


The  partial  derivative  terms  determine  the  impact  of  a  change  in 
the  control  on  the  two  adjacent  control  intervals.  However,  Eqs.  (4.29)  and 
(4.30)  are  both  valid  only  for  interior  node  points.  For  boundary  nodes,  these 
equations  break  down  since  at  either  boundary,  u<_i  or  0,4.1  does  not  exist. 
Therefore,  over  the  first  interval  only  Eq.  (4.30)  is  valid  while  over  the  last 
interval  only  Eq.  (4.29)  is  used.  Elsewhere,  the  partial  derivative  is  zero. 


With  the  differential  equations  for  $  and  \P,  and  the  general  gain 
equation  (Eq.  3.23),  the  gain  matrices  for  the  single  and  multiple  optimization 
trajectory  are  determined.  Because  the  matrix  has  zeroes  in  the  first 
column,  the  gain  matrix,  K,  also  has  zeroes  in  the  first  column.  These  zero 
gains  are  not  stored  in  either  the  single  or  multiple  trajectory  optimization. 


Intuitively,  the  zero  value  of  these  gains  makes  sense.  The  gains  in 
this  column  relate  deviations  in  t  to  control  changes.  As  the  final  value  of  t 
is  not  constrained,  its  variation  from  the  optimal  trajectory  does  not  influence 
the  control.  Conversely,  a  deviation  from  the  nominal  trajectory  for  y  and/or 
v  will  have  an  impact  on  the  desired  final  conditions  so  a  control  change  is 
necessary. 


For  both  methods  eight  gain  matrices  are  computed  from  u  =  0  up  to, 
and  including,  ti  =  0.875.  At  ti  =  0.875  there  are  two  nodes  and  two  boundary 
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constraints.  At  the  last  node  (u  =  1)  there  are  not  enough  nodes  to  satisfy 
the  boundary  conditions.  Also,  since  the  last  node  is  at  the  final  constraint 
manifold,  it  would  take  an  infinite  control  to  correct  for  perturbations. 

4.3  Simulation 

To  test  both  gain  calculation  methods,  simulations  of  the  lunar  launch 
problem  are  made.  Each  simulation  uses  a  fourth-order  Runge-Kutta  integrator 
to  simultaneously  integrate  the  differential  equations  of  motion  for  a  perturbed 
and  unperturbed  trajectory.  The  simulation  integrates  the  equations  of  motion 
using  40  steps  for  a  step  size  of  Au  =  0.025.  The  perturbed  trajectory  has  either 
lunar  gravity  or  thrust  acceleration  perturbed  ±5%  from  their  nominal  values. 
The  perturbed  and  unperturbed  trajectories  are  sampled  at  each  integration 
step  to  determine  the  error  vector,  6x.  Simulations  are  done  closed-loop  using 
a  sample  and  hold  technique  for  the  control  update.  On  each  integration  step 
the  nominal  control  vector  is  linearly  interpolated  while  the  control  update 
is  held  constant.  To  provide  a  basis  for  comparison,  another  simulation  using 
open-loop  control  is  also  made.  The  open- loop  method  uses  the  optimal  control 
in  the  presence  of  the  perturbation  with  no  correction  to  the  nominal  control 


vector. 


Chapter  5 


Results 

5.1  Single  Trajectory  Optimization  and  Simulation 

For  the  single  trajectory  optimization  the  nonlinear  programming 
code  VF02AD  is  used  to  calculate  the  nine  node  suboptimal  control  starting 
at  u  =  0  up  to,  and  including,  8  =  1.  VF02AD  returns  the  following  results: 


Mi  I 

'  0.4539471496  ' 

26.00925579  ' 

02 

0.4068478813 

23.31066650 

03 

0.3580472572 

20.51459670 

<?4 

0.3079740487 

17.64561319 

05 

0.2592985250 

rad.  = 

14.85671112 

06 

0.2077640929 

11.90400566 

07 

0.1567470292 

8.98094322 

08 

0.1048252859 

6.00604647 

.  09  . 

.  0.0528867033  . 

.  3.03018489  . 

deg. 


(5.1) 


With  this  control  history,  the  vehicle  reaches  the  following  final  conditions  at 
Uf  =  5444  ft/sec  on  the  unperturbed,  nominal  trajectory: 


tf  =  272.7061  sec  (5.2) 

yj  =  50000.0005  ft  (5.3) 

vf  =  1.3987xl0-6  ft/sec.  (5.4) 


The  gain  program  uses  the  suboptimal  control  history  and  the  gain 
equation  to  compute  the  optimal  control  gains.  The  gain  vectors  to  three 
significant  figures  are  shown  in  Table  5.1. 
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NODE 

u 

NODE 

y 

GAINS 

t> 

1 

0.000 

-0.369E-5 

-0.673E-3 

2 

0.125 

-0.289E-5 

-0.462E-3 

3 

0.250 

-0.385E-5 

-0.521E-3 

|  4 

0.375 

-0.573E-5 

-0.640E-2 

j  5 

0.500 

-0.940E-5 

-0.831E-2 

6 

0.625 

-0.179E-4 

-0.118E-2 

7 

0.750 

-0.461E-4 

-0.201E-2 

8 

0.875 

-0.267E-3 

-0.581E-2 

Table  5.1:  Single  trajectory  optimal  control  gains  for  nodes  1  through  8 


To  check  the  accuracy  of  the  gains,  a  numerical  derivative  method  is 
used.  The  original  gain  matrix  at  each  node  relates  the  error  in  the  state  (dif¬ 
ference  between  the  perturbed  and  nominal  trajectory)  with  the  change  in  the 
control  at  each  node  along  the  trajectory,  required  to  keep  the  perturbed  tra¬ 
jectory  in  the  neighborhood  of  the  extremal.  This  relationship  can  be  expressed 
as 


6a  =  K6x.  (5.5) 

As  a  result,  K  can  be  approximated  as 

K  »  (5.6) 

Table  5.2  shows  a  comparison  of  the  gain  matrix  at  node  1  (before  eliminating 
extraneous  gains)  computed  with  the  gain  equation  and  numerical  partials. 

Most  of  the  differences  in  the  values  of  the  numerical  partials  can 
be  attributed  to  accuracy  and  round  off  errors  in  VF02AD.  To  compute  the 
derivatives,  VF02AD  calculates  the  optimal  control  based  on  a  perturbation  in 
the  state  vector.  The  state  pertrubation  is  10-4  while  the  convergence  accuracy 
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|  GAIN  MATRIX  FOR  ti  = 

_ « 

fj  GAIN 

NUMERICAL 

1  EQUATION  1 

PARTIALS 

-0.369E-5 

-0.673E-3 

-0.353E-5 

-0.676E-3 

-0..M3E-5 

-0.559E-3 

-0.283E-5 

-0.558E-3 

-0.192E-5 

-0.439E-3 

-0.193E-5 

-0.439E-3 

-0.980E-6 

-0.313E-3 

-0.999E-6 

-0.313E-3 

-0.265E-8 

-0.183E-3 

-0.117E-6 

-0.185E-3 

0.100E-5 

-0.480E-4 

0.110E-5 

-0.446E-4 

0.203E-5 

0.898E-4 

0.200E-5 

0.904E-4 

0.306E-5 

0.230E-3 

0.252E-5 

0.228E-3 

0.410E-5 

0.370E-3 

0.511E-5 

0.369E-3 

Table  5.2:  Numerical  partial  approximation  of  control  gains  for  node  1 


for  VF02AD  is  10-9.  This  produces  significant  figures  on  the  order  of  10-5. 
Therefore,  for  extremely  small  gains  (i.e.  -.265E-8)  the  numerical  partials  do 
not  agree  favorably  with  the  calculated  result. 

Tables  5.3  and  5.4  show  the  simulation  results  using  the  single  opti¬ 
mization  gains.  Each  table  lists  the  error  in  the  final  values  for  y  and  v  from 
the  desired  final  conditions.  To  provide  an  additional  basis  for  comparison,  the 
open-loop  results  have  also  been  included. 

Simulations  using  single  optimization  gains  show  a  marked  improve¬ 
ment  in  the  final  conditions  over  the  open-loop  results.  Positional  errors  due 
to  acceleration  perturbations  are  better  handled  than  errors  in  gravity,  how¬ 
ever.  Velocity  errors  at  u/  are  approximately  the  same  regardless  of  whether 
the  perturbation  is  due  to  gravity  or  thrust  accelerations. 

Figs.  5.1  and  5.2  depict  the  history  of  the  position  and  velocity  error 
as  a  function  of  the  variable  of  integration,  ti.  When  viewing  the  graphs  it 


%  Change  State  Deviation  from  Optimal 
in  a  Closed  Loop  Open  Loop 


10.530 


19.760 


21.84  +5.0 


8328.003 


-76.358 


-17.592 

6590.498 

+0.811 

2.987 

Table  5.3:  Single  optimization  method  for  ±5%  modeling  error  in  thrust  accel 
eration 
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Figure  5.1:  Single  optimization,  federation  perturbation 
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(i&) 

%  Change 
in  g 

State 

Deviation  from  Optimal  | 

Closed  Loop 

Open  Loop 

5.054 

l 

Cn 

O 

y 

65.178 

9891.024 

V 

-2.705 

72.540 

1  5.586 

+5.0 

y 

-63.989 

-9891.023 

L__ 

V 

2.616 

-72.540 

Table  5.4:  Single  optimization  method  for  ±5%  modeling  error  in  lunar  gravity 


Figure  5.2:  Single  optimization,  gravity  perturbation 
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must  be  remembered  that  the  objective  of  the  guidance  law  is  to  operate  in 
the  neighborhood  of  the  extremal  path.  As  a  result  the  position  and  velocity 
errors  will  not  normally  converge  to  zero  until  later  in  the  trajectory. 
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5.2  Multiple  Trajectory  Optimization  and  Simulation 

The  multiple  trajectory  optimization  method  also  uses  a  nine-node 
optimal  trajectory  as  the  basis  for  the  nominal  control.  A  nine  node  optimiza¬ 
tion  is  then  computed  from  each  node  on  the  suboptimal  trajectory  to  the  final 
constraint  manifold.  Table  5.5  shows  the  resulting  optimal  control  vectors  from 
each  node.  As  expected,  the  optimal  control  from  node  1  is  the  same  as  the 
single  trajectory  optimization. 

Table  5.6  lists  the  optimal  control  gains  at  each  node.  Compared  with 
the  single  optimization  gains,  the  multiple  optimization  are  up  to  twice  as  large 
as  gains  produced  under  the  single  optimization  method.  Gains  at  either  ends 
of  the  trajectory  are  approximately  equal. 

To  verify  the  accuracy  of  the  gains,  a  numerical  partial  approximation 
is  made  similar  to  that  discussed  in  the  previous  section.  Table  5.7  shows  that 
the  numerical  partials  agree  extremely  well  the  gains  computed  using  the  gain 
equation. 

Tables  5.8  and  5.9  and  Figs.  5.3  and  5.4  present  the  results  of  the 
simulation  for  the  multiple  optimization  case. 

On  the  whole,  comparisons  of  the  multiple  and  single  optimization 
methods  produced  mixed  results.  For  acceleration  perturbations,  single  opti¬ 
mization  gains  gave  slightly  better  results  than  multiple  optimization  methods. 
On  the  other  hand,  end  point  conditions  for  the  multiple  optimization  approach 
in  the  face  of  gravity  perturbations  were,  relatively  speaking,  much  better  than 
single  optimization  results. 

Figs.  5.3  and  5.4  depict  the  history  of  the  error  of  y  and  v  from 
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NODE  1 

0  =  0 

NODE  2 

0  =  0.125 

NODE  3 

0  =  0.250 

NODE  4 

0  =  0.375 

0.45394715 

0.40634061 

0.35766840 

0.30858725 

0.40684788 

0.32136083 

0.27792256 

0.35804726 

0.32106646 

0.28359996 

0.24597137 

0.30797405 

0.27670197 

0.24516890 

0.21360186 

0.25929852 

0.23386572 

0.20826123 

0.18270982 

0.20776409 

e mmsM 

isircr;T?3:n 

0.15027212 

0.15674703 

0.14389432 

0.13098104 

0.11795366 

0.10482529 

0.09830906 

0.09183870 

0.08527328 

0.05288670 

wmmnm 

0.05300320 

NODE  5 

0  =  0.500 

NODE  6 
u  =  0.625 

NODE  7 

0  =  0.750 

NODE  8 
ti  =  0.875 

0.25871156 

0.20780274 

0.15648354 

0.10482639 

0.23386278 

0.18893294 

0.14341997 

0.09789380 

0.20798410 

0.16950645 

0.13168537 

0.09302267 

0.18175316 

0.14960138 

0.11570712 

0.08265920 

0.15684631 

0.13073685 

0.08325916 

0.13100243 

0.11166393 

0.08838203 

0.06894193 

wiiwism.™ 

0.09213030 

0.08020577 

0.06716968 

0.07866887 

0.07212200 

0.06518209 

0.05876025 

0.05283867 

0.05343331 

0.05345789 

Table  5.5:  Optimal  control  vectors  for  nodes  1  through  8  (radians) 


NODE 

u 

NODE  GAINS  p 

y 

t; 

|  1 

0.000 

-0.369E-5 

-0.673E-3 

2 

0.125 

-0.494E-5 

-0.780E-3 

1  3 

0.250 

-0.688E-5 

-0.921E-3 

4 

0.375 

-0.101E-4 

-0.112E-2 

5 

0.500 

-0.161E-4 

-0.141E-2 

6 

0.625 

-0.290E-4 

-0.190E-2 

7 

0.750 

-0.661E-4 

-0.288E-2 

8 

0.875 

-0.267E-3 

-0.580E-2 

Table  5.6:  Multiple  optimization  control  gains 
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GAIN 

EQUATION 

NUMERICAL 

PARTIALS 

-0.369E-5 

-0.673E-3 

-.353E-5 

-.C75E-3 

-0.494E-5 

-0.780E-3 

-.495E-5 

-.781E-3 

-0.688E-5 

-0.921  E-3 

-.688E-5 

-.922E-3 

-0.101E-4 

-0.112E-2 

-.101E-4 

-.112E-2 

-0.161E-4 

-0.141E-2 

-.160E-4 

-.141E-2 

-0.290E-4 

-0.190E-2 

-.290E-4 

-.190E-2 

-0.661E-4 

-0.288E-2 

-.656E-4 

-.286E-2 

-0.267E-3 

-0.580E-2 

-.267E-3 

-.577E-2 

Table  5.7:  Numerical  partial  approximation  of  control  gains 

nominal  as  a  function  of  u  for  perturbations  due  to  thrust  acceleration  and 
lunar  gravity.  Both  figures  show  that  the  multiple  optimization  gains  provide 
good  control  of  the  error  and  an  almost  asymptotic  approach  toward  zero  error 
at  uj. 

The  thrust  acceleration  perturbation  graph  (Fig.  5.3)  illustrates  an 
apparent  inconsistency  due  to  the  change  of  the  variable  of  integration.  In 
the  graphs,  an  error  in  velocity  does  not  change  sign  while  the  position  error 
oscillates  once  around  the  zero  error  axis.  The  reason  for  this  has  to  do  with 
normalizing  the  equations  of  motion  with  u/  instead  of  tj. 

Changing  the  variable  of  integration  to  u  means  the  position  vector 
becomes  a  function  of  u  instead  of  t.  Therefore,  to  compute  an  error  vector, 
the  state  vector  of  the  perturbed  trajectory  at  some  up  is  compared  with  the 
state  vector  of  the  nominal  trajectory  at  the  same  up.  These  two  points  along 
the  trajectory,  however,  do  not  occur  at  the  same  point  in  time.  As  a  result, 
while  an  element  of  the  perturbed  state  vector,  for  example  yp,  may  be  greater 
than  the  nominal  trajectory  at  up,  yp  may  still  be  lower  than  it  should  be  for 
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(*) 

%  Change 
in  a 

State 

Deviation  from  Optimal  | 

Closed  Loop 

Open  Loop  | 

|  19.760 

-5.0 

y 

13.757 

-8328.003  I 

1 

V 

-3.630 

-76.358  I 

|  21.84 

+5.0 

y 

-27.812 

6590.498  | 

1 

V 

4.300 

69.086  1 

Table  5.8:  Multiple  optimization  method  for  ±5%  modeling  error  in  thrust 
acceleration. 


Figure  5.3:  Multiple  optimization,  thrust  acceleration  perturbation. 


32 


(.■£?) 

%  Change 
in  g 

State 

Deviation  from  Optimal  | 

Closed  Loop 

Open  Loop 

1  5.054 

-5.0 

y 

47.558 

9891.024 

1 

V 

-1.767 

72.540 

5.586 

+5.0 

y 

-47.295 

-9891.023 

V 

1.741 

-72.540 

Table  5.9:  Multiple  optimization  method  for  ±5%  modeling  error  in  lunar 
gravity. 
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Figure  5.4:  Multiple  optimization,  gravity  perturbation. 
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Figure  5.5:  Altitude  versus  time 

the  same  point  in  time.  Fig.  5.5  shows  the  plot  of  the  y  element  as  a  function  of 
time  for  a  -5%  error  in  thrust  acceleration.  As  can  be  seen,  the  -5%  perturbed 
trajectory  does  not  exceed  the  nominal  trajectory. 


Chapter  6 


Conclusions 


The  purpose  of  this  research  is  to  develop  a  more  optimal  method 
of  computing  the  control  gains  at  each  node  through  multiple  optimization. 
Multiple  optimization  reoptimizes  the  trajectory  from  each  control  node  on  the 
suboptimal  trajectory.  Each  subtrajectory  from  a  node  to  the  final  constraint 
manifold  has  the  same  number  of  control  nodes,  r,  as  the  original  suboptimal 
trajectory  and  represents  the  optimal  r-node  solution.  The  suboptimal  control 
vector  for  each  subtrajectory  is  used  to  compute  the  control  gains.  These  gains 
are  more  aligned  with  the  most  optimal  path  from  a  node  to  the  final  constraint 
manifold. 

To  test  the  control  law,  simulations  are  made  using  a  representative 
optimization  problem,  the  lunar  launch  problem.  The  lunar  launch  problem 
involves  placing  a  vehicle  into  lunar  orbit  in  minimum  time  subject  to  con¬ 
straints  on  the  final  states.  The  simulation  integrates  the  equations  of  motion 
for  both  the  perturbed  and  unperturbed  (nominal)  trajectory.  The  value  of 
each  state  vector  at  each  integration  step  is  sampled  and  used  to  determine  the 
error  vector.  Multiplying  the  error  vector  by  the  control  gains  gives  the  control 
update  to  be  applied  to  the  perturbed  trajectory. 

Neighboring  extremal  guidance  using  gains  calculated  by  multiple  op¬ 
timization  methods  produce  results  comparable  to  the  single  trajectory  opti- 
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mization  method.  For  perturbations  of  ±5%  in  either  gravity  or  thrust  acceler¬ 
ation,  endpoint  conditions  using  multiple  optimization  gains  never  exceed  more 
than  48  feet  in  position  and  4.3  ft /sec  in  velocity  from  the  desired  final  con¬ 
ditions.  Single  optimization  gains  give  nearly  the  same  results  with  endpoint 
conditions  never  varying  more  than  65  feet  and  3  ft/sec  from  the  desired  final 
conditions.  Comparisons  of  each  method  for  a  specific  model  perturbations  give 
mixed  results,  with  multiple  optimization  delivering  smaller  errors  in  the  final 
conditions  for  perturbations  due  to  gravity  and  single  trajectory  optimization 
performing  better  for  thrust  acceleration  errors. 

These  observations  are  made  cognizant  that  the  sample  problem  in 
this  research  does  not  define  an  acceptable  error  for  the  final  conditions.  Pre¬ 
vious  research  [1]  has  used  the  criterion  that  the  percentage  error  in  the  final 
states  should  not  exceed  the  percentage  error  in  the  perturbation.  Unfortu¬ 
nately,  for  errors  in  y/  this  criterion  seems  rather  liberal  and  is  easily  met  by 
both  the  single  and  multiple  optimization  method.  Conversely,  a  5%  error  in 
vj  is  zero,  a  little  too  restrictive. 

One  area  for  further  research  is  to  explore  the  impact  of  increasing 
the  number  of  control  nodes  along  the  trajectory.  Increasing  the  number  of 
control  nodes  and,  hence,  the  number  of  gains  decreases  the  the  gap  between 
the  point  at  which  the  last  gain  matrix  exists,  and  the  final  constraint  manifold. 
This  approach  should  improve  guidance  near  the  end  of  the  trajectory  and 
subsequently  decrease  the  error  at  the  final  constraint  manifold. 


Bibliography 


[1]  Hull,  D.  G.  and  Helfrich,  C.  E.,  “Neighboring  Extremal  Guidance  For 
Systems  With  A  Piecewise  Linear  Control.”  AIAA  Third  Aerospace  Planes 
Conference,  Orlando,  FI.  Paper  AIAA-91-5067 

[2]  Bryson,  A.  E.  and  Ho,  Y.  C.,  Applied  Optimal  Control.  Hemisphere  Pub¬ 
lishing  Corp.,  Washington  D.  C.,  1975. 

[3]  Atkinson,  L.  V.  et  al.,  Numerical  Methods  with  FORTRAN  77.  Addison- 
Wesley  Publishing  Co.,  Wokingham,  England,  1989 

[4]  Hull,  D.  G., “Optimal  Control  Theory”,  unpublished  notes,  Department  of 
Aerospace  Engineering  and  Engineering  Mechanics,  University  of  Texas, 
Austin,  Texas 


36 


VITA 


Michael  James  Nowak  was  born  August  2, 1957  in  Hampton,  Virginia, 
the  son  of  Donald  and  Mary  Nowak.  He  graduated  from  Rome  Catholic  High 
School  in  1975.  Mr.  Nowak  graduated  with  a  degree  in  Astronautical  Engi¬ 
neering  and  received  his  commission  as  a  second  lieutenant  from  the  United 
States  Air  Force  Academy  in  1979.  Following  graduation  Mr.  Nowak  spent 
five  months  as  a  research  assistant  at  the  Air  Force  Academy  working  on  the 
Global  Positioning  Satellite.  For  the  next  eleven  years  he  was  an  Air  Force 
pilot  and  instructor  in  training  and  tactical  reconnaissance  aircraft.  He  began 
his  masters  program  at  the  University  of  Texas  in  June  1991. 


Permanent  address:  409  RIDGEWOOD  DR. 

ROME  NEW  YORK  13440 


This  thesis  was  typeset1  with  DTgX  by  the  author. 


i 

OT]gX  document  preparation  system  was  developed  by  Leslie  Lamport  as  a  special  version 
of  Donald  Knuth’s  TfeX  program  for  computer  typesetting.  TfejX  is  a  trademark  of  the 
American  Mathematical  Society.  The  UTjrjX  macro  package  for  The  University  of  Texas  at 
Austin  thesis  format  was  written  by  Khe-Sing  The. 


